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SUMMARY 


This  report  discusses  the  work  accomplished  during  the  second 
quarter  (March  -  May  1963)  on  "Research  on  Unsteady  Stall  of  Axial 
Flow  Compressors",  Contract  No.  DA.  49-186-AMC-13(X),  for  the 
Harry  Diamond  Laboratories,  Washington,  D.  C.  ,  by  Cornell  Aeronautical 
Laboratory,  Inc. 

Initial  effort  in  this  program  is  concentrated  on  continuation  of 
previous  research  at  CAL  concerning  (1)  an  experimental  study  of  the 
separation  of  a  laminar  boundary  layer  from  walls  moving  upstream  and 
downstream  in  a  steady  flow  using  a  shrouded  rotating  cylinder  in  a  wind 
tunnel,  and  (2)  an  analytical  study  of  compatibility  conditions  between 
circulation  disturbances  on  an  actuator  disk  in  a  mean  flow  with  swirl  and 
the  appearance  of  vorticity  downstream  of  the  disk. 

Hot  wire  anemometer  equipment  previously  on  order  was  received. 
This  equipment  was  modified  to  provide  proper  temperature  compensation 
for  use  with  nonstandard  probes.  Boundary  layer  velocity  profiles  were 
obtained  over  approximately  half  the  range  of  angular  positions  on  the 
cylinder  which  appear  to  be  necessary  in  order  to  define  separation  velocity 
profiles  for  upstream-  and  downstream-moving  walls. 

The  analysis  for  compatibility  Conditions  governing  the  existence 
of  rotating  circulation  disturbances  on  a  single-stage  axial-flow  compressor 
disk  was  continued.  The  very  severe  restrictions  imposed  on  the  results 
previously  reported  (radial  perturbation  velocity  component  assumed  equal 
to  zero,  no  radial  variation  in  flow  variables)  were  relaxed  to  the  extent 
that  only  the  assumption  that  the  radial  perturbation  velocity  component 
was  assumed  zero.  Under  this  assumption  and  imposing  the  conditions 
that  the  disturbance  flow  be  bounded  with  time  and  axial  distance  from  the 
compressor  disk,  it  was  found  that,  as  before,  the  only  possible  rotating 
circulation  disturbance  on  the  compressor  disk  is  one  rotating  with  the 
blades. 
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II  INTRODUCTION 


Rotating  stall  has  been  a  continuing  problem  for  the  designers  of 
axial-flow  turbojet  engines  since  the  inception  of  their  use  in  aircraft.  For 
the  compressor  designer,  the  major  problem,  associated  with  rotating  stall 
has  been  one  of  blade  fatigue  failure  resulting  from  blade  resonance  with 
the  passage  of  stall  zones.  From  a  performance  point  of  view,  unsteady 
stall  phenomena  may  place  severe  limitations  on  the  acceleration  or  decelera¬ 
tion  of  a  jet  engine. 

Extensive  research  effort  has  been  devoted  to  this  problem  during 
the  past  two  decades.  Progress  was  made  in  that  a  degree  of  understanding 
of  the  phenomenon  was  attained;  however,  the  basic  problem  remained 
unsolved  when  lack  of  funds  reflecting  a  declining  interest  in  jet  engine 
development  brought  this  research  effort  nearly  to  a  halt  in  the  late  1950's. 

In  particular,  the  unsteady  boundary -layer  separation  phenomena  which, 
presumably,  underlie  rotating  stall  of  compressor  elements  continue  to  pose 
problems  of  fundamental  interest. 

Due  to  a  revival  of  interest  in  axial-flow  jet  engine  development,  this 
program  was  initiated  to  continue  a  methodical  approach  to  solutions  for 
the  problems  posed  by  unsteady  compressor  stall.  The  work  reported  here 
is  a  continuation  of  a  research  program  concerning  rotating  stall  in  axial- 
flow  compressors,  carried  out  by  the  Cornell  Aeronautical  Laboratory  over 
a  three-year  period  ending  in  January  1959.  The  ultimate  objective  of  this 
research  effort  is  an  increased  understanding  of  rotating  stall  and  other 
unsteady  surge  characteristics  of  axial  compressors.  Included  will  be 
studies  of  the  aerodynamics  of  compressor  stages,  studies  of  unsteady 
sectional-lift  characteristics  of  single-blade  elements  near  stall,  and  studies 
of  possible  methods  of  alleviating  unsteady  surge  characteristics.  In  the 
current  year's  effort,  the  emphasis  is  on  the  continuation  of  two  studies 
started  under  the  original  program:  (1)  a  theoretical  study  of  the  role  of 
entering  swirl  pattern  in  causing  rotating  stall  of  a  single  stage,  and  (2)  an 
experimental  study  of  the  separation  criteria  for  a  moving  surface  (specifi¬ 
cally  a  shrouded  rotating  cylinder). 
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WORK  ACCOMPLISHED 


A.  Rotating  Cylinder  Experiments 

Details  of  the  motivation  of  this  work  and  of  the  basic  experimental 
apparatus  being  employed  were  reported  in  the  First  Quarterly  Progress 
Report  (Reference  1)  and  will  not  be  repeated  here. 


Hot  Wire  Anemometer  Calibration 


The  Flow  Corporation  hot  wire  air  velocity  meter,  Model  55B1, 
ordered  during  the  first  quarter  (purchased  with  CAL  funds)  was  delivered 
early  in  the  period.  This  unit  was  purchased  with  an  adapter  (Model  A.-3) 
which  permitted  the  use  of  a  standard  Flow  Corporation  HWP-A  probe  in 
single-wire  operation  with  the  meter  in  place  of  the  special  two-wire  probe 
supplied  with  the  Model  55B1  unit  (HWP-A  probes  were  available  at  CAL 
from  the  previous  rotating  stall  research  program).  This  adapter  was 
modified  to  permit  use  of  two  HWP-A  probes  in  two-wire  operation  (for 
temperature  compensation)  after  preliminary  calibration  of  the  single-wire 
unit  showed  excessive  scatter.  After  some  experimentation  with  methods 
for  mounting  wires  on  the  probes  to  obtain  good  solder  joints,  consistently 
good  calibration  data  were  obtained.  Figure  1  shows  calibration  curves  for 
the  first  three  hot-wires  used  during  the  boundary-layer  velocity  surveys. 

Preliminary  boundary -layer  velocity  profile  data  obtained  with  an 
HWP-A  probe  indicated  considerably  thinner  boundary-layer  on  the  rotating 
cylinder  than  the  data  previously  obtained  with  a  total-head  probe.  The 
total-head  probe  data  is  probably  not  reliable  for  the  inner  portions  of  the 
boundary  layer  near  the  cylinder  surface  where  displacement  effects  and 
large  shear  gradients  would  tend  to  distort  the  velocity  profile.  Such  data 
should  be  reliable,  however,  near  the  outer  edge  of  the  boundary  layer. 

It  was  concluded  that  the  HWP-A  probes  were  causing  flow  interference  in 
the  boundary  layer.  Therefore,  a  special  boundary-layer  probe  was  pur¬ 
chased  (Flow  Corporation  Model  HWP-D).  Velocity  profiles  obtained  with 
this  probe  compared  well  with  total-head  probe  data  for  the  outer  portion 
of  the  boundary  layer.  All  boundary-layer  data  reported  here  have  been 
obtained  with  this  boundary-layer  probe. 


Boundary- Layer  Velocity  Profiles 


Figure  2  shows  the  notation  appropriate  to  the  boundary-layer 
velocities  (note  that  this  notation  is  different  from  the  notation  shown  in 
Figure  1(b)  of  Reference  1). 

=  the  free-stream  velocity. 

U  -  the  outer-flow  velocity  on  the  edge  of  the  boundary  layer. 
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Figure  I  HOT  WIRE  PROBE  CALIBRATIONS 
(TNO-WIRE  MODE  OF  OPERATION) 
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u(y)  =  boundary  layer  velocity  component  parallel  to  the  local 
tangent  to  the  cylinder  surface. 

2/YyJ  =  boundary  layer  velocity  component  perpendicular  to  the 
local  tangent  to  the  cylinder  surface. 


Velocity  profiles  in  the  boundary  layer  on  the  rotating  cylinder 
have  been  obtained  in  the  range  of  O  from  78°  to  17.0°  in  increments  of 
10°  for  SI  =  +2000,  +1500,  tlOOO,  +500,  +200  and  0  rpm. 

The  hot-wire  probe  is  affected  by  immediate  proximity  to  the 
surface  of  the  cylinder  (within  0.  010  inch).  This  effect  is  apparently 
associated  with  heat  transfer  between  the  hot  wire  and  the  cylinder  surface. 

It  was  found  possible  to  correct  the  experimental  data  to  account  for  this 
wall  proximity  effect,  and  the  correction  appears  to  be  negligible  for  most 
of  the  data,  becoming  of  mildly  significant  magnitude  only  very  close  to 
the  wall  at  Q  -  120°  and  fl  =  0  and  -200  rpm. 

A.  further  difficulty  is  posed  by  the  fact  that  the  hot  wire  measures 
the  magnitude  of  the  velocity  vector  normal  to  the  wire;  hence,  in  the 
region  of  the  boundary  layer  where  velocities  parallel  to  the  cylinder  surface 
approach  zero,  the  (very  small)  velocity  component  perpendicular  to  the 
cylinder  surface,  Wy)  ,  is  recorded  by  the  hot  wire.  It  is  in  this  area 
that  the  data  are  of  most  interest  and  importance  as  regards  the  velocity 
profiles  at  separation.  There  appears  to  be  no  way  to  experimentally  deter¬ 
mine  the  magnitude  of  this  component  with  the  existing  equipment,  nor  is 
there  available  theory  for  moving-wall  boundary  layers  of  use  in  estimating 
the  magnitude  of  the  vertical  velocity  component  near  separation.  To  over¬ 
come  this  difficulty,  various  assumptions  on  the  shape  of  the  V" (y)  profile 
have  been  used  to  correct  the  experimental  data.  The  results  are  illustrated 
in  Figure  3  for  6  =  100°  and  A  =  -513  rpm.  It  is  apparent  that  the 

corrected  profile  is  relatively  insensitive  to  the  shape  of  the  assumed  pro¬ 
file  for  v(y)  .  This  is  generally  true  of  all  the  boundary-layer  data 
obtained  to  date.  Thus,  it  can  be  hoped  that  the  corrected  profiles  are 
good  approximations  to  the  true  profile  of  the  tangential  velocity  component 
in  the  boundary  layer. 

Velocity  profiles  for  0,  +500  and  -500  rpm  are  shown  in  Figures 
4,  5  and  6.  In  Figure  4,  the  0  rpm  velocity  profiles  behave  as  expected 
with  an  apparent  separation  point  {diL/dy^.0  ~  0)  at  9  555  120°.  Velocity 
profiles  for  0  >  120°  have  not  as  yet  been  obtained.  The  boundary-layer 
profiles  in  Figure  5  are  very  thin  (note  the  expanded  scale  for  y  ),  and 
there  is  no  indication  of  separation  being  approached  at  the  maximum  value 
of  &  .  The  profiles  shown  in  Figures  4  and  5  contain  no  correction  for 

WyJ  • 


The  profiles  obtained  at  -500  rpm  (Figure  6)  have  been  corrected 
for  the  perpendicular  velocity  component  by  assuming  that  'V'(y)  =  constant 

and  is  equal  in  magnitude  to  the  minimum  velocity  recorded  by  the  hot  wire, 
The  boundary  layer  in  this  case  grows  much  faster  than  the  0  rpm  case,  but 
there  is  no  indication  of  an  abrupt  increase  in  thickness  which  would  indicate 
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Figure  5  BOUNDARY-LAYER  VELOCITY  PROFILES  AND  THICKNESSES  vs 


Figure  6  BOUNDARY-LAYER  VELOCITY  PROFILES  AND  THICKNESSES  vs  6 


separation  by  the  accepted  standard  for  zero  wall  velocity.  The  profiles 
at  large  values  of  6  appear  to  be  flattening  against  the  y  -axis  as 
hypothesized  in  Reference  2,  but  the  abrupt  change  from  a  predicted  rela¬ 
tively  large  flattened  portion  near  the  wall  to  the  wall  velocity  at  the  wall 
is  not  apparent.  The  profiles  appear  to  approach  the  wall  velocity  in  a 
fairly  gradual  fashion  with  decreasing  distance  from  the  wall.  The  behavior 
of  the  profiles  shown  here  and  the  gradual  increase  in  the  boundary-layer 
thickness  implied  by  the  smoke  pictures  of  Reference  3  suggest  that  separa¬ 
tion  does  not  occur  at  a  well-defined  point  but  is,  rather,  a  quite  gradual 
process.  Experimental  data  to  be  obtained  at  values  of  9  greater  than 
120°  should  help  to  clarify  this  idea. 

B.  Theoretical  Results 


The  initial  theoretical  effort  during  this  program  is  concerned  with 
a  continuation  of  the  work  reported  in  Reference  4.  This  work  was  a  study 
of  possible  disturbance  patterns  which  may  exist  in  the  flow  field  produced 
by  a  single-eompr essor  stage  --with  a  view  to  the  description  of  possible 
three-dimensional  or  partial-span  modes. 

It  should  be  made  clear  that  this  theoretical  study  is  not  directed 
toward  attempts  to  predict  specific  characteristics'  of  rotating  stall  phenomena, 
such  as  speed  of  propagation,  number  of  stall  zones,  or  the  prediction  of 
when  rotating  stall  will  occur  in  terms  of  compressor  characteristics. 

Rather,  the  investigation  is  concerned  with  compatibility  conditions  in  the 
flow  which  determine  whether  rotating  stall  is  a  permissible  disturbance  of 
the  particular  type  of  flow  considered. 

Moore  {Reference  4)  showed  that  for  small-disturbance  axial 
flow  with  mean  swirl  of  the  free-vortex  type  through  an  axial  compressor 
disk,  the  only  type  of  rotating  circulation  disturbance  on  the  compressor 
disk  that  could  exist  is  one  independent  of  blade  radius  which  is  fixed  and 
rotating  with  the  blades;  for  this  type  of  disturbance,  the  resulting  unsteady 
flow  downstream  of  the  compressor  was  irrotational.  The  theoretical  work 
of  the  present  program  is  an  attempt  to  extend  Moore's  analysis  to  more 
general  types  of  swirl  in  the  steady  flow  upstream  of  the  compressor  disk. 

It  is  shown  that  if  the  radial  component  of  disturbance  velocity  is  zero  (or 
more  generally  if  the  derivative  of  this  component  with  respect  to  axial 
distance  is  zero)  for  any  mean-flow  swirl  the  only  rotating  circulation  dis¬ 
turbance  that  can  exist  on  the  blades  must  rotate  with  the  blades  if  the 
resulting  disturbance  flow  is  to  remain  bounded  with  both  increasing  time 
and  with  axial  distance  downstream  of  the  compressor.  (Moore's  results 
apply  whether  the  disturbance  flow  remains  bounded  or  not.  )  Under  these 
conditions  the  resultant  unsteady  flow  is  irrotational. 

In  particular,  it  would  appear  that  these  results  would  be  applicable 
to  flow  through  an  axial  compressor  with  a  narrow  annulus  where  it  is  reason¬ 
able  to  expect  that  radial  disturbance  velocities  would  be  very  small. 

Linearized  equations  for  the  convection  of  vorticity  were  developed 
in  References  4  and  1,  for  a  cylindrical  coordinate  system  with  origin  at 
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the  center  of  the  compressor  disk:  *  is  the  axial  coordinate,  positive 
downstream;  r  is  the  radial  coordinate,  and  6  is  the  angular  coordinate. 
These  equations  are 


(J-+U  A  +  UoSA){ 

\H  U°dx  +  r  i>&lj 


dx 


-mr*>riv+-u?(rs)'— 


r  Lr  ‘rsx 


Asr5)^+,-(rS)'--s]ri 


(1) 


where  U0  is  the  free  stream  mean  axial  velocity, 

U0S  is  the  tangential  mean  velocity, 

(<5,  0ATe  the  perturbation  components  of  vorticity,  and 
(U,V'I4U) are  the  perturbation  velocity  components. 


i 


The  perturbation  vorticity  components  are  related  to  the  perturbation 
velocity  components  by 

v  -  J-  _  J_  dtr 

’  ”  p  dr  r  30 

J_  c?u,  _  t 

?  =  r  39  ax 

r  -  _  A u 

’  ~  d%  dr 

Also,  the  continuity  equation  is 

du,  1_  a(rv')  J_  a^or  _ 
a  x,  r  a  r  r  a& 


(2) 


(3) 


These  equations  are  subject  to  the  requirement  of  periodicity  with  respect 
to  9  in  the  velocity  and  vorticity  components. 

If  it  is  assumed  that  VJ  =  O  ,  Equations  (1)  become 
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ra(rS)  ^  r  d(rS) 
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where 


X  -  J-  CMdll 

h  ~  r  dr 

-  1  du  B-iM- 

(  r  de  dK 


c  = 


-  da 
dr 


and  Equation  (3)  becomes 


da_  J_  SsUs  _ 
dx  rd&~ 


A.  solution  for  ^  is  then 

-  g(\ ,r) 


/ 


(4) 


(5) 


(6) 


(7) 


where 

A  -  9  -  (a,  +  ^)x +- &J0t 

and  ^U0  is  the  rotational  velocity  of  the  circulation  disturbance  on  the 
compressor  disk.  Note  that  periodicity  in  9  for  n  requires  that 
be  periodic  in  A 


Consider  next 


i£ 

d-t 


+  u  £& 

0  dx 


Do 5  d£  Ug  $d(r£)  d*a/  .  f B(rS) 

r  d9  r  l  dr  dx  [dr 


(8) 


from  Equations  (4).  In  particular,  consider  the  particular  integral  which 
satisfies  the  nonhomogeneous  term  (A/',  £ d/$r (r S ) - 5 ]  g  (  A  ,  r  )  .  Assume  this 

particular  integral  has  the  form  U0/r[ci(%,  \- ,  0 ,  t  j  g  fA  ,r) 


Substituting  in  the  left-hand  side  of  Equation  (8),  and  equating  with 
U°Jr[d/dr(rS) -S  ]  <j(X,  r  )  there  results 


a  =  ae(X,r)  + 


'd(rs) 

dr 


-5 


( 1-ou,  U0)i+CL < *] 


The  function  Qa  (X,r)g(X,r)  can  be  incorporated  into  the  complementary 
solution  of  Equation  (8).  In  Equation  (9),  there  is  only  one  arbitrary  con¬ 
stant;  hence,  it  can  be  required  that  f  be  bounded  as  X,  becomes  large, 


(9) 
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or  as  ~t  becomes  large,  but  not  both.  Hence,  if  f  is  to  be  bounded 
for  both  jc  and  -d  ,  it  must  be  required  that  q(\tr)  ~  n  -  0 


From  Reference  4,  Equation  (45),  the  disturbance  circulation  / 
is  related  to  the  appearance  of  radial  vorticity  at  the  compressor  disk  by 


^  -t 

3t  '  ° 


i/  s  U 


(10) 


where  ~k(J0  is  the  angular  velocity  of  the  blades  and  ^ 

Hence,  if  q  =  0,  then 

=  H(e-itu0t)  +  F(r)  (11) 

only.  Therefore,  for  bounded  small  perturbation  flows  with  swirl  a  general 
function  of  radius  in  the  mean  flow,  the  only  rotating  circulation  disturbance 
possible  is  one  fixed  in  the  blades,  if  VJ  =  0;  note  that  this  is  true  if  it  is 
only  required  that  dlr/3x,  -  0,  i.  e.  ,  that  ^  is  not  a  function  of  % 

At  the  compressor  disk, 

=  r-Jl-hS2^  =  F'(r) 


as 

3  °> 

(Reference  4,  Equations  (46)  and  (38).  )  Hence,  the  unsteady  small  disturbance 
flow  downstream  of  the  compressor  is  irrotational,  as  the  flow  upstream 
is  irrotational,  and  the  unsteady  circulation  disturbance  results  in  no  vorticity 
shed  into  the  flow  downstream. 

It  is,  of  course,  possible  that  even  if  the  requirement  of  bounded¬ 
ness  for  the  small  disturbance  flow  had  been  relaxed,  the  same  conclusion 
would  have  been  reached,  as  was  the  case  for  the  free-vortex  swirl  flow 
treated  in  Reference  4;  this  case  has  not  as  yet  been  treated  in  the  current 
analysis. 
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IV 


WORK  FOR  THE  NEXT  QUARTER 


A..  Experimental  Program 

It  is  expected  that  the  rotating  cylinder  experimental  work  will  be 
concluded  during  the  next  quarter.  A  heavy  workload  is  scheduled  for  the 
CA.L  One -Foot,  High-Speed  Wind  Tunnel  in  which  this  work  is  being  per¬ 
formed.  It  is  desirable,  therefore,  to  complete  our  tests  as  soon  as 
possible  to  avoid  excessive  delays  and  cost  associated  with  disassembly 
and  reassembly  of  the  apparatus  that  might  result  from  fitting  our  tests 
into  the  tunnel  schedule. 


B.  Theoretical  Program 

Work  will  continue  on  the  compatibility  analysis  currently  underway. 
Effort  will  initially  be  directed  toward  extending  the  analysis  by  removing 
the  restriction  that  2 ir/ dx  be  assumed  zero. 


C.  General 

Consideration  will  be  given  to  the  most  desirable  areas  for  future 
work  on  the  program. 
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